MET receptor serves as a promising target in melanoma brain metastases

The development of brain metastases hallmarks disease progression in 20–40% of melanoma patients and is a serious obstacle to therapy. Understanding the processes involved in the development and maintenance of melanoma brain metastases (MBM) is critical for the discovery of novel therapeutic strategies. Here, we generated transcriptome and methylome profiles of MBM showing high or low abundance of infiltrated Iba1high tumor-associated microglia and macrophages (TAMs). Our survey identified potential prognostic markers of favorable disease course and response to immune checkpoint inhibitor (ICi) therapy, among them APBB1IP and the interferon-responsive gene ITGB7. In MBM with high ITGB7/APBB1IP levels, the accumulation of TAMs correlated significantly with the immune score. Signature-based deconvolution of MBM via single sample GSEA revealed enrichment of interferon-response and immune signatures and revealed inflammation, stress and MET receptor signaling. MET receptor phosphorylation/activation maybe elicited by inflammatory processes in brain metastatic melanoma cells via stroma cell-released HGF. We found phospho-METY1234/1235 in a subset of MBM and observed a marked response of brain metastasis-derived cell lines (BMCs) that lacked druggable BRAF mutations or developed resistance to BRAF inhibitors (BRAFi) in vivo to MET inhibitors PHA-665752 and ARQ197 (tivantinib). In summary, the activation of MET receptor in brain colonizing melanoma cells by stromal cell-released HGF may promote tumor self-maintenance and expansion and might counteract ICi therapy. Therefore, therapeutic targeting of MET possibly serves as a promising strategy to control intracranial progressive disease and improve patient survival. Supplementary Information The online version contains supplementary material available at 10.1007/s00401-024-02694-1.


Introduction
The interaction of brain colonizing tumor cells with the tumor microenvironment (TME), mainly comprising innate and adaptive immune cells, microglia, astrocytes, neurons and oligodendrocytes crucially determines the developmental stages of brain metastases (BM).Brain metastases are observed in 20-40% of melanoma patients during the course of disease and micrometastases are evident in more than 75% of autopsied brains [51].Hence, only a subset of melanoma cells that entered the brain develop symptomatic and detectable BM during the lifetime of melanoma patients.Unlike peripheral metastases, the emergence of BM depends on a plethora of environmental cues such as the spatiotemporal availability of factors that are provided by cells of the TME, supporting or repressing tumor cell growth [65].Moreover, single-cell RNA sequencing (scRNAseq) studies have confirmed regional heterogeneity of astrocytes [29], oligodendrocytes [27] and microglia [36,67] in healthy human brains.Particularly, astrocytes and microglia adopt a reactive cell state [6,38] that accompanies the secretion of pro-and anti-inflammatory factors under pathological conditions [42,67].It is therefore possible that subfractions of astrocytes and microglia communicate and react with tumor cells in different ways.Probably, neuroinflammation precedes colonization of the brain by tumor cells.However, tumor cells invading the brain amplify inflammatory processes mediated by astrocytes and infiltrating tumor-associated microglia and macrophages (TAMs) [12,62].Recently, signaling mediated by hepatocyte growth factor (HGF) and the related receptor MET (c-MET, HGFR) was identified as the trigger of reactive microglia [52].HGF is thereby secreted by microglia in the context of trauma but also under normal conditions and seems to play a special role in the growth and self-renewal of neural stem cells in the subventricular zone (SVZ) of rat brains [45,73].Therefore, metastatic melanoma cells expressing MET might scavenge HGF from the brain for activation of processes downstream of MET mediating survival and proliferation.
Here, we used transcriptome and methylome profiling to unravel the epigenetic and transcriptomic landscapes of MBM that featured infiltration of TAMs with emphasis on the potential role of microglia in the activation of the HGF/MET receptor signaling pathway.The MET receptor inhibitors PHA-665752 and tivantinib (ARQ197) effectively blocked the growth of brain metastases-derived cells (BMCs).Hence, targeting of MET receptor signaling might serve as a potent therapeutic target for brain metastases lacking druggable BRAF V600 mutations.

Patient cohorts
All procedures performed in this study were in accordance with the ethical standards of the respective institutional research committees and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.All patients gave written informed consent for the collection and scientific use of tumor material which was collected at the Biobank of the Charité -Comprehensive Cancer Center (CCCC).The study was approved by the Ethics Committee of the Charité (EA1/152/10; EA1/107/17; EA4/028/18 and EA1/107/17 and EA1/075/19) and Universitätsmedizin Greifswald (BB 001/23).

Live cell imaging-based drug sensitivity assays
Drug treatments were performed 24 h after seeding of 2500-5000 cells/96-well in 100 µl medium.The response of BMCs and conventional melanoma cell lines to dabrafenib, PHA-665752 or ARQ197 (all purchased from Selleckchem) in a range of 1 nM-10 µM of eight technical replicates was determined by live cell imaging.Images were taken every three hours using a 10 × objective and the general label-free mode, two pictures of eight technical replicates per condition were taken.Drug response was assessed by changes in the cellular density over time.The cell density was determined by a confluence mask tool as part of the IncucyteS3 software.IC50 values were calculated by curve-fitting (https:// search.r-proje ct.org/ CRAN/ refma ns/ REAT/ html/ curve fit.html) based on confluence measurements on day 3.

In vivo experiments
All animal experiments were performed in accordance with the German Animal Protection Law under the permission number G0130/20 obtained via the Berlin Ministry of Health and Social Affairs (LaGeSo).ARRIVE 2.0 Guidelines were strictly followed and performed as previously reported [49].Briefly, 2.5 × 10 4 BMC1-M4 and BMC2 cells were stereotactically inoculated into brains of female Crl:CD1-Foxn1 nu nude mice (8-9 weeks of age, 24-26 g, Charles River Laboratories) were with using a 1 µl Hamilton syringe and a stereotactic frame as described previously [2].Tumor growth was tracked by MRI and animals were sacrificed by perfusion with 4% PFA in deep anesthesia after tumors reached a volume of 20 mm 3 .Following, whole brains were removed, dehydrated, paraffin embedded and sections of 2 µm were used for downstream analyses.

RNA isolation and sequencing
Isolation of total RNA from snap-frozen tumors and RNA sequencing was performed as previously reported [49].Briefly, 100 ng of total RNA was used for library preparation with TruSeq Stranded Total RNA Sample Preparation-Kit and Ribo-Zero Gold Kit (Illumina).Paired-end (2 × 100 bp) sequencing of RNA libraries with integrity numbers (RIN) ≥ 7 was performed on the NovaSeq6000 platform at Cegat GmbH, Tuebingen (Germany).Following demultiplexing of sequenced reads and adapter trimming [33], FASTQ files were obtained.Raw counts of proteincoding genes were normalized using the DESeq2 (https:// bioco nduct or.org/ packa ges/ release/bioc/html/DESeq2.html) package [40].Differential expression of genes between groups was determined after fitting models of negative binomial distributions to the raw counts.Raw p-values were FDR (false discovery rate)-adjusted for multiple testing and a value below 0.05 for the adjusted p-values were used to determine significant differentially expressed genes.

Quantitative real-time RT-PCR
RNA isolation from frozen cell pellets was performed with the RNeasy Mini Kit (Qiagen, Germany) and, following the manufacturer's protocol as previously reported [49].qRT-PCR was carried out on a Step one plus PCR cycler (Applied Biosystems, Germany) for 30-40 cycles.Primers were designed for 55-60 °C annealing temperatures.Relative expression levels were calculated with the ΔΔCT method [39], normalized to β-actin.Primer sequences are shown in Supplementary Table 7.

A microglia-specific gene cluster discriminates MBM
Microglia are a unique population of antigen-presenting cells in the central nervous system (CNS) that are capable of clearing the brain of microbes, dead cells and protein aggregates [59].Besides, microglia play a crucial role during injury repair and display an exceptional role in immune surveillance and tumor clearance [9,13].Although the role of tumor-associated microglia and macrophages (TAMs) in primary brain tumors such as glioblastoma [3,8,68,69] has been intensively studied, their role in the progression of brain metastases remains poorly understood.
We performed immunohistochemistry (IHC) of our MBM cohort (Supplementary Table 1 and [49]) to determine the levels of activated TAMs expressing Iba1 (AIF1).Although Iba1 serves as a well-established marker, reactive microglia cannot be distinguished from brain-infiltrated macrophages [34].We observed that Iba1 levels classified MBM into highly and lowly TAM infiltrated tumors (Fig. 1a, Supplementary Fig. 1a).Moreover, we observed overlapping patterns of infiltration of Iba1 high TAMs and CD3 + T cells (Fig. 1b).As CD3 only provided information about levels of T cell infiltration, we used the ESTIMATE algorithm [75] to gain insight into the overall degree of immune cell infiltration of MBM.In line with our previous observation, tumors with intensive TAM and T cell infiltration exhibited a high immune score (Pts 3, 4, 10, 12) whereas MBM with low levels of Iba1 high / CD3 + cell infiltration (Pts 1, 2) or low expression of Iba1 (Supplementary Fig. 1b, c) showed low immune scores (Fig. 1c).As expected, brain metastases derived cell lines (BMCs) with absence of immune cells featured lowest scores (Supplementary Fig. 1b, c).The brain has long been considered a sanctuary where tumor cells can grow undisturbed and protected from attack by immune cells.We observed a high correlation between Iba1 and the immune score by analyzing the expression levels of Iba1 in an independent dataset of brain (MBM) and extracranial metastases (EM) (Fig. 1d), suggesting a relationship between the degree of infiltration of TAMs and immune cells not only in the brain.
As high levels of immune (T) cell infiltration are generally associated with good prognosis [55], we determined the probability of survival related to Iba1 expression of patient´s with (study EGAS00001003672) and without (TCGA-SKCM) MBM.We observed beneficial effects of high Iba1 levels in the TCGA cohort (HR = 0.46 (0.35-0.62), logrank p = 1.3e-07) (Fig. 1e, f), but observed no beneficial effect on the survival of MBM patients.Since no data on TAMinfiltrated MBM are available, we performed comparative methylome and transcriptome profiling of a core set of Iba1 high (n = 5, methylome profiling; n = 10, transcriptome profiling) and Iba1 low (n = 2, methylome profiling; n = 6, transcriptome profiling) tumors which were selected based on IHC.We identified a set of 416 differentially methylated genomic regions (DMRs) that corresponded to 294 MBM expressed genes (Fig. 1g) a core set of markers (n = 31) sufficient to split tumors (Fig. 1h; Supplementary Table 2, 3).Among them, we identified the integrin family member   Tumor TILs and gut-homing receptor ITGB7 -which we described in our previous study as a marker distinguishing BRAF and NRAS mutant MBM [49]-and APBB1IP (amyloid b precursor protein-binding family b member 1 interacting protein).Both are associated with better prognosis in patients with colorectal cancer [19,80] and clustered with known TAM-associated genes such as P2RY12 and AIF1 (Fig. 1h).Remarkably, all clustered tumors were associated with a high immune score.A correlation analysis of clustered genes revealed a high degree of correlation among each other (Fig. 1i) and association with hepatocyte growth factor (HGF) that was recently connected with microglia activation [52].However, only some of the identified markers within the gene cluster such as APBB1IP were specifically expressed in microglia but not in brain-infiltrating macrophages or other brain cells (Fig. 1j).APBB1IP has been identified as a conserved microglial gene [20] and binding partner of amyloid precursor protein (APP), Tau, 14-3-3 g, and glycogen synthase kinase 3 b (GSK3 b) was associated with actin dynamics and retinoic acid signaling [31,35] and was significantly (MBM: R = 0.86, p < 2.2e-16) correlated with immune score (Fig. 1k) and survival of melanoma patients (Supplementary Fig. 1d-e).Moreover, our survey identified a differentially methylated side (Supplementary Table 4) within the promoter of PD-L2 (PDCD1LG2) that may predict progression-free survival in melanoma patients receiving anti-PD-1 immunotherapy [28].PD-L2 expression was associated with favorable survival (p = 0.020) of patients with MBM (Supplementary Fig. 1f).We found additional genes among our cluster that were expressed in TAMs and significantly associated with immune score (Supplementary Fig. 1 g-n).

Expression of ITGB7 serves as indicator of immune cell infiltration
Recent studies have shown that ITGB7 plays a critical role in the recruitment of T cells to the intestine and that downregulation of ITGB7 is important in protecting intestinal tumors from attack by activated T cells [10,80].Hence, we sought to investigate ITGB7 in more detail.Mining of publicly available immune cell data (studies GSE146771 [79], DICE database [61]) revealed expression of ITGB7 across different immune cell stages including naïve and memory subsets of T cells, B cells and NK cells (Fig. 2a and Supplementary Fig. 2a).We found that ITGB7 was rather expressed in MBM with infiltration of immune cells and particularly within immune cell dense areas (Supplementary Fig. 2b).Co-staining revealed accumulation of CD3 + T cells as well as of Iba1 high TAMs (Fig. 2b).In line, we found higher levels of CD4, CD274 (PD-L1) and Sushi Domain Containing 3 (SUSD3) in MBM featuring high ITGB7 expression (Fig. 2c) and validated a potential, previously observed [49] correlation of ITGB7 and SUSD3.ITGB7, SUSD3 and APBB1IP showed expression across different immune cell types except for monocytes and NK cells (Supplementary Fig. 2c-f).Global (850 k) methylome profiling uncovered four epigenetic regulation sites of ITGB7 (Supplementary Table 4) with two sites that were associated with gene expression and immune score (Fig. 2d, left and center panel, Supplementary Fig. 3a), located in a proximal enhancer-like region (probe cg26689077) or nearby the promotor of ITGB7 (probe cg01033299).The latter site was also identified in the TCGA-SKCM cohort.The sites did not correlate with the BRAF mutation status of MBM (Fig. 2d, right panel) in contrast to additional two sides that were found within intergenic regions including an CpG island located between exons 4 and 5 (probes cg11510999 and cg18320160; Supplementary Fig. 3b-e).We suggest that ITGB7 may serve as an indicator of the degree of immune A recent study demonstrated that MBM feature a lower T cell content than matched extracranial metastases, however, response rates to ICi of both were comparable [71].Assuming that ITGB7 expression might be crucial for T cell recruitment, we ascertained the levels in MBM (n = 79) and EM (n = 59; study EGAS00001003672).ITGB7 was expressed in both metastatic subtypes and was significantly correlated (MBM: R = 0.51, p = 1.8e-06;EM: R = 0.69, p = 1.1e-09) with the tumor´s immune scores (Fig. 2e).As we suggest that ITGB7 expression might indicate the degree of immune cell infiltration and possibly serve as indicator of response to ICi, we next performed correlation analysis of ITGB7 and known markers of T cells and B cells.We observed a high concordance with immune cell-related but not tumor cell-related genes (NGFR, MITF, MLANA, SLC45A2) and correlation with expression of PDCD1LG2 and SUSD3, irrespective of the side of metastasis (Fig. 2f,  g).In line with previous observations, ITGB7 was expressed in primary and metastatic tumors (TCGA-SKCM) and like SUSD3 was associated with favored survival (Fig. 2i).In summary, our survey identified a set of markers that are potentially associated with the level of TAM/immune cell infiltration, particularly ITGB7 might serve as a marker for a favorable course of the disease.

A signature-based deconvolution revealed MET receptor signaling in microglia-enriched MBM
Our previous survey identified a set of markers that potentially characterize a molecular subset of MBM, likely showing a favorable course and response to ICi therapy [23,30].To further characterize the molecular subsets, we performed a single-sample Gene set Enrichment-Analysis (ssGSEA) using signatures that defined immune-related signaling or processes that involved MET receptor or STAT3 signaling (Supplementary Table 5).We found that MBM with high immune score were enriched in genes associated with MET and STAT3 signaling, tumor inflammation, stress and senescence (SenMayo [58]) (Fig. 3a) and featured the presence of reactive microglia, astrocytes and immune cell subsets, among them stem cell-like CD8 + T cells (TCF7) [46] in tumors, absent in BMCs.CD8 + (TCF7) T cells are necessary for long-term maintenance of T cell responses and predicted positive clinical outcomes [14, [57].Signatures clearly discriminated MBM and BMCs and reinforced the differences between Iba1 high (Pts 3, 4) and Iba1 low/neg (Pts 1, 2) tumors.We therefore suggest that the activation of MET-or STAT3-mediated signaling processes or those related to stress/senescence or inflammation strongly depends on the composition of the tumor microenvironment, likely determining the response to therapeutic interventions.Although infiltration of TAMs is not evident in all MBM, microglia infiltration seems to be an early occurring process observed ~ 21d after intracranial injection of BMCs into brains of immune-compromised Crl:CD1-Foxn1 nu mice [49] (Fig. 3b).The activation of Stat3 signaling in tumoradjacent cells (Fig. 3c), may propose a rapid response of brain microenvironmental cells to brain colonizing tumor cells.We observed a comparable pattern of enrichment of molecular processes in Iba1 high MBM of an independent (study EGAS00001003672 [17], n = 79 MBM) (Supplementary Fig. 4a).
HGF or scatter factor (SF) is the only identified ligand of MET and exerts pivotal functions during neural development, regulating the growth and survival of neurons [15, Fig. 3 Signature-based deconvolution identified the parameter of MBM featuring a favorable disease course and identified a role of MET signaling.a Single-sample GSEA (ssGSEA)-based deconvolution of MBM of study EGAS00001005976 using customized gene signatures indicating "Signaling" processes, cellular subsets and stages of microglia and astrocyte and immune cell subsets.ssGSEA demonstrated distinct separation of MBM with high, median or low immune score regarding expression levels of signature genes, BMCs served as controls.ssGSEA uncovered differentially activated pathways and processes such as MET and STAT3 and interferon signaling, senescence (SenMayo), stress response and tumor inflammation in tumors enriched for reactive microglia and astrocytes and innate and acquired immune cells subsets.b Confocal microscopy images of orthotopic tumors established by stereotactic injection of BMC1-M4 or BMC2 cells into brains of Crl:CD1-Foxn1 nu mice [49], stained for Iba1 (green, microglia) or Iba1, GFAP (red, astrocytes) and KBA.62 (turquoise, pan-melanoma cell marker).DAPI served as a nuclear counterstain.Markers show distinct areas of tumor (MBM) and microenvironment (TME) and regions of microglia infiltration, 21 days after intracranial injection [49].MBM-TME boarders are indicated by white, dashed lines.c IHC of tumors investigated in (b) for activation and tyrosine phosphorylation (residue Y705) of STAT3.pSTAT3 Y705 is particularly present in microenvironmental cells (astrocytes).Black, dashed lines indicate MBM-TME boarders.In b, c, bars indicate 50 µm.d-e Expression levels of hepatocyte growth factor (HGF) in tumors of studies EGAS00001005976, TCGA-SKCM and EGAS00001003672 demonstrating HGF expression in all tumor subsets.f, g Investigation of HGF expression in immune cell subsets (DICE database [61]) and brain cells (study GSE73721) revealed the highest levels in basophil granulocytes and monocytes (f) and in astrocytes and microglia (g).h UMAP projection of expression profiles from nuclei isolated from 5 neurotypical donors as provided by Seattle Alzheimer's disease brain cell atlas (https:// portal.brain-map.org/ explo re/ seatt le-alzhe imers-disea se), cellular subtypes are color coded (left panel).Log-normalized expression levels of HGF in nuclei isolated from 5 neurotypical donors (center panel).Log-normalized expression levels of HGF in nuclei isolated from 84 aged donors (42 cognitively normal and 42 with dementia), right panel, demonstrating an increased number of HGF expressing microglia and astrocytes as triggered by inflammatory processes.i Dot plot showing the correlation of HGF expression and immune score of BM (R = 0.49, p = 5.3e-06) and EM (R = 0.41, p = 1.5e-03) indicating a potential role of HGF in immune cell-related processes.Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (d-g).Significance was determined by unpaired, two-sided t-test (e) or one-way ANOVA (g) , likely serving as an inducer of reactive microglia by an autocrine loop in response to trauma or neurodegenerative disorders [52].Therefore, MET-expressing melanoma cells infiltrating the brain can benefit from the HGF-regulated systems that naturally occur in the brain and employ them as a survival strategy.We observed HGF expression among tumors of different data sets comprising MBM, EM and primary tumors (studies EGAS00001005976; TCGA-SKCM; EGAS00001003672) with no significant difference in HGF levels of tumor subsets (Fig. 3d, e).Investigation of immune cell and brain cell data (DICE database [61] and study GSE73721 [78]) revealed high expression of HGF in monocytes and astrocytes (Fig. 3f, g), and the exploration of single-cell studies (GSE115978 [32] and GSE186344 [22]) revealed melanoma and MBM associated microglia/ macrophages as a source of released HGF (Supplementary Fig. 4b).
Assuming that the mutual interaction of tumor cells and TAMs determines the routes of MBM progression and fosters the activation of MET receptor-related processes, we investigated levels of MET signaling-associated genes.We found that degrees of HGF, PIK3CG, PTK2B, STAT3 and MAP4K1 significantly correlated with microglia score (Supplementary Fig. 4c-e) that was defined as average expression (log2 FPKM) or methylation status (β-value) of microglia markers APBB1IP, SYK, HCK and P2RY12 (Supplementary Table 6).HGF might be released by homeostatic and reactive microglia (RM) [1,52] or reactive astrocytes (RA) [37,64] but transcriptome profiling data on MBM-related RM/RA are not available.We surveyed the Seattle Alzheimer´s Disease Brain Atlas which is implemented in the Allen brain atlas database (https:// portal.brain-map.org/).Dementia fostered the expansion of microglia and astrocytes with increased expression of HGF (Fig. 3h, center and right panels).Reactive microglia and immune cell released HGF might hence be responsible for the activation of growth factor/survival signaling in adjacent tumor cells.The level of HGF expression significantly correlated with an immune score in the brain (BM, R = 0.49, p = 5.3e-06) and extracranial metastases (EM, R = 0.41, p = 1.5e-03), (Fig. 3i).

Expression and activation of MET receptor classifies a molecular subset of MBM
Understanding the molecular mechanisms that establish cellular dependencies and thus control the development and maintenance of brain metastases is critical for their therapeutic manipulation.Recently, we identified that the expression of Ecad and NGFR sufficiently discriminated molecular subsets of MBM [49] likely to exhibit different responses to therapeutics and mechanisms of interaction with cells in the microenvironment (Fig. 4a).To identify potential druggable targets, we surveyed the pan-MBM, NGFR and Ecadspecific gene sets for cell surface receptors that may serve as crucial key factors controlling tumor cell survival and maintenance.We identified 24 receptors that distinguished Ecad + and NGFR + tumors (Fig. 4b).Particularly ADIPOR1 (adiponectin receptor 1, p = 1.9e-02),SIRPA (signal regulatory protein alpha, p = 1.1e-05) and PLXNC1 (plexin C1) showed significantly increased expression in Ecad + subsets of MBM and EM (Supplementary Fig. 5a).In addition, we found MET receptor predominantly expressed in Ecad + tumors (p = 1.4e-04) and significantly (p = 2.7e-05) higher expressed in MBM than EM (Fig. 4c, left and center panels).MET was enriched in proliferating tumor cells featuring high levels of markers of cell cycle progression such as cyclins b1 and b2 (CCNB1, CCNB2), proliferating-cell nuclear antigen (PCNA) and Ki67 (MKI67) (Fig. 4c, right panel, Supplementary Table 6) and likely define yet another subset of MBM.
We investigated the distribution of MET expression and of Ecad (CDH1) and MITF, the transcriptional regulator of MET in melanocytes as well as marker of pigmented melanoma cells among primary tumors, extracranial metastases (TCGA-SKCM), and brain metastases (MBM_CHA).

◂
We observed distinct as well as shared tumor cell states (Fig. 4d).Moreover, we detected that the MET receptor was not only present in Ecad + /MITF ± tumors but rather showed gradual distribution across tumors and co-expression in a minority of NGFR + cells (Fig. 4d and Supplementary Fig. 5b).To additionally unravel the cell states, we next explored single-cell transcriptome data (study GSE115978) for distribution of MET expression in melanoma cell subsets comprising immune cells, cancer-associated fibroblasts (CAF), endothelial cells, macrophages, tumor cells and natural killer (NK) cells.We observed that MET is primarily expressed in tumor cells and to a lesser extent in immune (T cells) cells/NK cells, macrophages, CAF or endothelial cells (Supplementary Fig. 5c, d).Tumor cells showed discrete and overlapping expression of MET, Ecad and MITF (Supplementary Fig. 5d-f) and were co-expressed with NGFR in a small subset of tumor cells (Supplementary Fig. 5 g, h).The investigation of BMCs (BMC1-M1) validated the cooccurrence of MET and NGFR expression in cellular subsets (Supplementary Fig. 5i, j), suggesting that the MET tyrosine kinase receptor pathway may serve as a potent survival and maintenance mechanism in MBM and that the therapeutic targeting by small molecule inhibitors [77] may eliminate several tumor cell fractions including NGFR + cells potentially featuring resistance to BRAFi.
Next, we assessed whether expressed MET (Fig. 4e) indeed participated in active signaling processes.Phosphorylation of MET at tyrosine residues 1234/1235 (pMET Y1234/1235 ) is critical for kinase activation and initiation of downstream processes and was evident in nearly all MET high MBM investigated, independent of the BRAF mutation status (Fig. 4f, g).MET receptor alterations are evident in 9% of all SKCM melanoma cases, including amplification as observed in 1.13 -17.2% or 11% of melanoma (TCGA-SKCM, study by Ramani et al. [50]).However, targeted DNA sequencing (TargetSeq) and fluorescence in-situ hybridization (FISH) revealed the absence of MET activating mutations and a tendential MET amplification in only one case (Pat 5, Supplementary Fig. 6a, b).However, all but one tumor (Pat 14) showed high polysomy.As we assume that microglia/macrophages may foster activation of MET receptor signaling in a subset of tumor cells, we performed co-IHC for pMET Y1234/1235 and Iba1.We observed pMET Y1234/1235 positive tumor cells in close proximity to Iba1 high TAMs (Fig. 4h), though the MET receptor was not activated in Iba1 high microglia that resided in adjacent normal tissue (Supplementary Fig. 6c, upper panel).However, MET receptor activation was also evident in scattered tumor cells in the absence of adjacent Iba1 high TAMs (Supplementary Fig. 6c, lower panel) suggesting paracrine mechanisms or additional sources of HGF such as immune cells or astrocytes.Considering that HGF levels, like those of other growth factors provided by stromal cells, might depend on spatial factors, we examined the Allen Brain Atlas database and found that HGF is comparably expressed in different brain sections (frontal lobe (FL), parietal lobe (PL), temporal lobe (TL), occipital lobe (OL)) but is lowly abundant in the brainstem (pons) (Supplemental Fig. 6d).The spatially dependent expression of growth factors in the brain may therefore determine the dependencies of the tumor cells.

Interferon signaling determines response of MBM to immune checkpoint inhibitor therapy
Interferon-gamma signaling has been identified as an important mechanism for the upregulation of PD-L1 on melanoma cells and escape from immune recognition.On the other hand, recent studies uncovered that high interferon-gammarelated gene expression signature scores (IFN-γ score) were associated with low risk of melanoma relapse from neoadjuvant ipilimumab plus nivolumab therapy [53,56].

The targeting of MET receptor serves as a promising strategy to control MBM growth
Although a fraction of MBM exhibits immune cell subset enrichment and interferon response signatures and responds to ICi therapy, MET-expressing brain metastatic melanoma cells may benefit from HGF released by stromal cells to drive progression.Hence, activation of MET signaling may depend on the degree of tumor-stroma interaction, possibly counteracting the beneficial impact of immune checkpoint inhibition (ICi).Resistance-mediating processes include the phosphorylation of ribosomal protein S6 (pS6), which is downstream of MET and mTOR signaling [16] and was observed in progressive BRAFi-resistant melanomas [63,74].
We assessed pS6 phosphorylation of serine residues 235/236 and found co-occurrence of activated MET receptor and of pS6 235/236 in MITF positive tumors (Fig. 5a and Supplementary Fig. 8a, b).Moreover, pS6 phosphorylation was evident in a BRAF wt (T2002) and mutated (V600E, BMC53) cell lines probably suggesting a general activation of pS6 signaling irrespective of the presence of mutated BRAF (Fig. 5b).As MET signaling might serve as mediator of cell survival and therapy resistance, we assessed the efficacies of the ATP-competitive inhibitor PHA-665752 (PHA) and the non-ATP-competitive, clinical phase I/II MET receptor inhibitor (METi) ARQ197 (tivantinib) in BMCs that showed variable levels of MET expression (Fig. 5c and Supplementary Fig. 8c).ARQ197 failed to improve the outcome and overall survival of patients with hepatocellular carcinoma [81] but may potentially be effective in melanoma patients.The initial testing revealed a general response of BMCs (BMC1-M1, BMC53), T2002 cells and conventional cell lines (A375, A2058, MeWo) to both inhibitors irrespective of the BRAF mutation status (Fig. 5d, e).Next, we investigated the response of BMCs that harbored druggable BRAF V600 mutations or non-druggable BRAF mutations (N581Y, BMC2) or cell lines with wildtype BRAF alleles (MeWo) to low (0.3-10 nM) and high doses (> 10 nM-10 µM) of PHA.We observed a general response of all cell lines to PHA (Fig. 5f) irrespective of their sensitivity to dabrafenib (Fig. 5g) or the general presence or absence of BRAF mutations.We observed IC 50 values within a range of 1.4 µM (A375 cells) to 2.5 µM (MeWo cells), validating the previous results.Considering the higher sensitivity of cell lines to ARQ197, we next evaluated the response of intrinsically resistant cell lines (BMC4, druggable BRAF V600 but refractory; BMC2) showing only moderate or no response to dabrafenib (Fig. 5g) as indicated by IC 50 values (BMC4, IC 50 = 226.4nM and BMC2, IC 50 = 3029.5nM) to low (0.3-10 nM) and high doses (> 10 nM-10 µM) of ARQ197.As for PHA, the broad range (1 nM-10 µM) testing of ARQ197 in BMCs, T2002 and conventional melanoma cells achieved a general response in all cell lines tested (Figure h-j, Supplementary Fig. 8d).However, we observed that the non-(brain) metastatic cell lines A375, A2058, and T2002 were more sensitive to ARQ197 treatment.The median IC 50 value of BMCs was ~ 600 nM (range: 406.5-800.1 nM).Cell lines lacking BRAF and NRAS mutations (MeWo, T2002) showed the highest responses to ARQ197 (Fig. 5k).Neither ARQ197 nor PHA have been clinically tested in patients with melanoma and brain metastases.The novel ATP-competitive and blood-brain-barrier-permeable inhibitor capmatinib is approved by the FDA for the treatment of NSCLC patients with hyperactivation of the MET pathway caused by an exon14 skipping mutation in the MET gene [69].We therefore tested the response of capmatinib in BMC, expecting a higher sensitivity as observed for PHA and ARQ197.However, unlike for NSCLC cell lines showing high sensitivity and against our expectations, BMCs neither responded to low (< 10 nM) nor high (> 1-10 µM) doses of capmatinib (Supplementary Fig. 8f) which is in line with the results of a recent study, suggesting a general low response of melanoma cell lines [4].
In summary, brain metastatic as conventional melanoma cell lines responded to METi, suggesting that targeting MET signaling might be a promising tool for the treatment of non-BRAF V600 and BRAF V600 mutated MBM that acquired resistance to BRAFi or for combinatorial of METi and ICi in NRAS mutated tumors and all MET expressing cellular subsets including BRAFi refractory NGFR + cells.

Discussion
The spatiotemporal development of primary and secondary brain tumors is strongly determined by the crosstalk of tumor and brain micronenvironmental cells, particularly macrophages, astrocytes and microglia [48] and the consequential activation of inflammatory processes [54].Although the neuro-inflammatory processes that are activated alongside the development and progression of primary brain tumors such as glioblastoma have been intensively studied, the mechanisms that accompany the emergence of brain metastases on the other side are not well investigated.
Here, we used combined transcriptome and methylome profiling to unravel the molecular features of MBM of different progression stages showing high and low level of tumor-associated macrophages/microglia (TAMs) infiltration, irrespective of the phenotype (Ecad, NGFR).Generally, TAMs foster the development and progression of primary brain tumors [25,76], however, their functional role in MBM may be different.We observed that MBM containing a high proportion of TAMs were associated with a high immune score and infiltration of CD3 + T cells.The profiling of Iba1/AIF1 high tumors revealed a cluster of genes, among them ITGB7, APBB1IP as SUSD3 and PD-L2, that were widely expressed among immune cell subtypes and previously associated with increased immune T cell infiltration [10,19,49,80] and favored outcome.Previous mouse studies demonstrated a pivotal role of Itgb7 for intestinal T cell recruitment and correlated low levels of Itgb7 with colorectal cancer progression and maintenance of intestinal stem cells via Ecad-mediated interaction [10,80].However, we observed clear protein expression of ITGB7 in tumors cells at immune cell-enriched areas of Ecad + and NGFR + tumors, suggesting a broader function in different subtypes of metastases and cancers.Enhanced expression of ITGB7 might be a prerequisite for immune cell invasion.Therefore, epigenetic marks that correlate with the expression of ITGB7 and other genes mentioned above may be of prognostic importance, and the expression of these markers could determine the pathways of intracranial progression.
As previously described for the Ecad + and NGFR + subtypes of MBM, whether tumors are enriched or depleted in TAMs and immune cell subsets is critical and may determine the response to therapeutic interventions.The subsequent ssGSEA-based characterization of MBM of studies performed by us and others revealed molecular programs fostering or accompanying the TAM + /TIL + tumor subtype.TAM + /TIL + tumors featured activation of MET and STAT3 signaling, increased stress response, tumor inflammation, senescence and activation of microglia and astrocytes but also activated interferon signaling.STAT3 activation in tumor-adjacent astrocytes in response to brain damage or tumor cells is well-investigated [26,64] and was rapidly induced in response to intracranially injected BMCs.Hence, the enrichment of STAT3 signature genes was likely attributed to tumor-adjacent astrocytes and infiltrated immune cells [47].However, STAT3 activation in brain metastases might foster brain colonization of melanoma cells [72].
The HGF/MET receptor signaling plays a pivotal role during brain development and neuro-regeneration, homeostasis of microglia and neurons [15,73] but is also involved in microglia activation in response to trauma [15,41,52].Brain infiltrating melanoma cells hence may engage the HGF/MET signaling axis of brain cells and utilize it for regulation of survival and proliferation.We observed the expression of MET receptor in the subset of E-cadherin (Ecad) expressing tumors [49], suggesting that Ecad + but not NGFR + cells may depend on HGF/MET signaling.Nevertheless, PCA representation of MET levels in individual tumors and single-cell resolution revealed a gradual expression pattern showing distinct MET high and MET low expressing subsets and overlapping expression with NGFR in tumor cells with intermediate levels of MET.This suggests that not only NGFR and Ecad but additional correlating genes may underlay cellular plasticity-related regulation of expression, hence therapeutic targeting of the MET receptor signaling pathway by small molecule inhibitors may eliminate several MET + tumor cell subsets.Immunohistochemistry indicated MET activation in MBM, occurring either in an autocrine or paracrine manner.Although we observed that MBM may provide HGF, whether and to which extent brain infiltrating melanoma cells provide HGF is unknown and needs additional investigation.As we observed that HGF is expressed by immune cell subsets and homeostatic or reactive astrocytes and microglia, we investigated the level of activated/ phosphorylated MET receptor in TAM-adjacent tumor cells.We found that tumor cells but not Iba1 high TAMs that ◂ resided in tumor cell-free adjacent stroma showed activation of MET, however, MET was also activated in the absence of adjacent microglia in some tumor cells, suggesting a paracrine effect of HGF.In line with our previous study [49], we observed the enrichment of interferon-response signatures in the subset of TIL high /immune score (IS) high tumors.We observed enrichment of interferon-response genes in MBM with high levels of ITGB7 expression and observed significant response of Itgb7 and Hgf among known interferon-inducible genes such as Cd274 [18] and Mx1 [69] in interferon-gamma treated BV2 murine microglia cells (unpublished study GSE132739).Hence, T cell-provided interferon-gamma might not only induce expression of CD274/PD-L1 but may also activate expression of HGF and ITGB7.Therefore, autocrine MET receptor signaling might be triggered in response to immune cell-released interferongamma and/or paracrine activation of MET signaling may occur via (INFG-activated) reactive glia-released HGF.
Our study bridges the gap between the immune cell phenotype of MBM and the activation of potentially therapeutic counteracting signaling pathways.The infiltration of TAMs and immune cells thus represents a double-sided sword and on the one hand is associated with an effective response to immune checkpoint inhibitors, but on the other hand can support the growth of MET-expressing tumor cells via secreted factors such as HGF.
Therefore, we finally assessed the potential role of small molecule inhibitors of MET receptor (METi) for targeting MBM that lack druggable BRAF V600 mutations or developed refractory disease.To this end, we took advantage of our well-characterized BMCs serving as in vitro model systems.We observed that the ATP-competitive inhibitor PHA-665752 [11] and the non-ATP-competitive, clinical phase II inhibitor ARQ197 (tivantinib) [44] elicited response in BMCs and conventional melanoma cell lines irrespective of the BRAF/NRAS mutation status.However, although being effective at doses of 100-200 nM in MeWo and A375 cells, ARQ197 showed a median IC 50 value of ~ 1 µM in BMCs, suggesting a general difference among brain metastatic and long-term maintained conventional cell lines established from either non-metastatic (A375) or locally metastatic (MeWo) tumors.Moreover, cell lines responded even less sensitive to PHA-665752 in a range of 1-2.5 µM and barely responded to capmatinib.This novel FDAapproved inhibitor features blood-brain barrier penetrance and hence might be a very promising drug for the treatment of brain metastases.However, being very effective in nonsmall cell lung cancer (NSCLC) patients and cell lines with MET amplification, marked MET overexpression or MET exon 14 skipping mutations and HGF-mediated activation of MET signaling, the low response of melanoma cell lines to capmatinib is in line with results of a recent study [4].The efficacy of capmatinib was recently tested in two clinical trials NCT03484923 and NCT02587650.However, the mechanisms of action of inhibitors specifically targeting tumors with hyperactivated MET e.g.caused by MET exon 14 skipping mutations are not fully elucidated.Using targeted sequencing, we were unable to detect MET exon 14 skipping mutations in our MBM cohort and MBM-derived cell lines in our previous study [49].Mechanistically, capmatinib and PHA-665752 act as ATP-competitive inhibitors, while ARQ197 blocks the activation of MET in an allosteric manner.The mechanistic difference in conjunction with the prevalent genetic aberrations such as activating BRAF mutations likely determines the response of tumor cells to MET inhibition.Although the inhibition of MET in MBM by brain-penetrable small molecule inhibitors is promising, the circumstances of MET activation in MBM and therapeutic interventions to selectively inhibit the MET receptor signaling pathway need to be intensively investigated in further studies.As ARQ197 was reported to block MET independently from MET binding [5], this suggests additional effects of the inhibitor.Hence we cannot exclude that ARQ197 caused cytotoxic effects that were not primarily associated with levels of MET expression or activation.
In summary, we have shown that MET receptor signaling is active in a subset of MBM, conferring a survival/ growth benefit independent of BRAF/NRAS mutation status.MET activation may occur in response to HGF released by TAM/immune cells and could counteract therapeutic interventions.Furthermore, we suggest interferoninduced expression of HGF in tumor cells triggered by interferon-gamma provided by stromal cells mediates autocrine activation of MET-signaling in tumor cells (Fig. 6).In addition, we demonstrated that methylome profiling of MBM has high potential to uncover gene regulatory sites such as identified in ITGB7, APBB1IP, SUSD3 and PD-L2 (PDCD1LG2) that may predict favorable progression of intracranial disease.

Limitations of the study
The present study is not without limitations.Due to the risks of a craniotomy and the partially limited regional accessibility of the tumors, alternative strategies such as radiotherapy or systemic therapy are the method of choice for patients with MBM.Therefore, only a small proportion of tumors are surgically removed, so that access to the corresponding tissue is considerably restricted.Although we suggested that HGF/MET signaling is activated in tumor cells in close proximity to infiltrated microglia, the efficacy of METi in established brain metastases/tumors remains to be proven particularly whether and to which extent METi affect normal homeostatic processes e.g.those crucial for neuron survival.

Fig. 1
Fig. 1 Transcriptome and methylome profiling of Iba1 high and Iba1 neg MBM revealed the identification of subset-specific genes.a Immunohistochemistry (IHC) for Iba1 (red) of MBM of indicated patients.b Representative IHC for levels of CD3 in Iba1 high (Pat 4) and Iba1 low/neg (Pat 1) MBM.c Immune score of MBM (study EGAS00001005976, n = 16) indicating different immunologic (color coded) subsets of tumors.d Dot plot showing the significant correlation of Iba1/AIF1 expression and immune score of brain metastases (BM, R = 0.86, p < 2.2e-16) and extracranial metastases (EM, R = 0.78, p = 5.5e-13).e Survival analysis of patients with MBM (study, EGAS00001003672), featuring high or low level of Iba1/AIF1 expression revealed no significant difference (p = 0.11).f Survival analysis of TCGA melanoma patients (n = 459), featuring a high or low level of Iba1/AIF1 expression revealed a significant difference (logrank p = 1.3e-07) and Cox-regression analysis showed association with favorable disease course (HR = 0.46).g Schematic representation of candidate identification by methylome and transcriptome profiling of n = 16 MBM of study.Methylome (850 k) profiling of Iba1 high (n = 5) or Iba1 low/neg (n = 2) identified 416 differentially methylated regions (DMRs), within the 5´-UTR of 316 corresponding genes of which 294 were expressed in MBM with 56 genes (77 DMRs), significantly (p ≤ 0.05) discriminating Iba1 high and Iba1 low/neg MBM.h Heat map representation of 77 DMRs (left panel) and top expressed (right panel) genes (n = 31).Analysis identified a panel of 12 genes that clustered with expression of microglia/TAM-associated genes AIF1, SYK and HCK.i Correlation analysis of cluster genes with association to immune/TAM regulated processes, the strength of the correlation is color coded.j Comparative t-SNE representation of brain cell subclasses microglia, neurons and oligodendrocytes (left) and expression of APBB1IP (Amyloid Beta Precursor Protein Binding Family B Member 1 Interacting Protein), expression level (log2 RPKM) is color coded.k Dot plot showing the significant correlation of APBB1IP expression and immune score of brain metastases (BM, R = 0.86, p < 2.2e-16) and extracranial metastases (EM, R = 0.92, p < 2.2e-16).Significance was determined by unpaired, two-sided t-test (d, g, k)

Fig. 2
Fig. 2 Expression of ITGB7 serves as an indicator of lymphocyte infiltration.a Box plot representation of levels of ITGB7 indicates a wide pattern of expression among indicated immune cell populations.Monocytes and neutrophil granulocytes show low levels of ITGB7.b IHC of a representative MBM of a patient with refractory intracranial disease for Iba1 (red, first column) and CD3 (brown, second column) indicating focal enrichment of microglia/macrophages and CD3 + T cells within ITGB7 positive areas (red, second column).Hematoxylin and eosin (H&E) staining shows discrimination of tumor cells and tumor-infiltrating lymphocytes (TILs) c Expression (FPKM, log2) of CD4, PD-L1 (CD274) and SUSD3 in MBM with a high or low level of ITGB7, indicating cellular co-occurrence.d Dot plot showing the significant inverse correlation (R = -0.87,p = 5.2e-05) of β-values (probe cg26689077) indicating the methylation level at a side located within the proximal enhancer-like structure of the ITGB7 gene and immune score of MBM (n = 14) of study EGAS00001005976 (first panel).Box plots represent a significant (p = 4.5e-04) or non-significant (p = 0.86) association of ITGB7 methylation (probe cg26689077) or BRAF mutation status (center and right panels) of all MBM investigated (n = 21).e Dot plot showing the significant correlation of ITGB7 expression and immune score of MBM (R = 0.51, p = 1.8e-06) and EM (R = 0.61, p = 1.1e-09) indicating immune-related expression of ITGB7 irrespective of the side of metastasis.f.) Correlation map showing high association (p < 0.05) of ITGB7 with relevant immune cell markers such as PD-1 (PDCD1), PD-L1 (CD274), PD-L2 (PDCD1LG2) but low correlation with tumor cell markers NGFR, MITF, MLANA or SLC45A2.g Dot plot showing the significant correlation of ITGB7 and expression of PD-L2 (BM: R = 0.45, p = 3.4e-05; EM: R = 0.42, p = 1.1e-03) and SUSD3 (BM: R = 0.44, p = 5.2e-05; EM: R = 0.61, p = 2.6e-07).h Dot plot showing the correlation of ITGB7 expression and immune scores of primary (PT; R = 0.59, p = 9.4e-16), metastatic (EM; R = 0.78, p = 2.2e-16) and brain metastatic (BM; R = 0.2, p = 0.61) melanoma (TCGA-SKCM), indicating that expression of ITGB7 is independent of melanoma progression stages.i Survival analysis of TCGA melanoma patients (n = 459), featuring a high or low level of ITGB7 and SUSD3 expression revealed a significant difference (log rank p = 4.0e-04 and p = 6.6e-08) and Cox-regression analysis showed association with favorable disease course (HR = 0.60 and HR = 0.48).Box and whisker plots show the median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (a, c, d).Significance was determined by unpaired, two-sided t-test (c, d) or one-way ANOVA (a) ◂ signaling Target of interferon signaling I mmune cells I nterferon-induced/ Innate iimmunity Interferon-induced/ immune cell response

Fig. 4
Fig.4 Ecad + MBM are defined by expression of MET receptor.a Schematic summary of the initial screen of MBM expression data of our recent study (EGAS00001005976; n = 16 MBM) for subset expressed receptors.MBM contains Ecad + and NGFR + subsets and admixed cells such reactive microglia, labeled by expression of Iba1/AIF1 and or P2RY12.The initial survey yielded 24 receptors that potentially establish cell survival/growth of MBM.b.) Correlation map (Spearman, p < 0.05) showing the relationship of identified receptors expressed in MBM of our previous study, emphasizing the distinct pattern of Ecad + and NGFR + molecular subsets.The value of the correlation coefficient is color coded.c Box plots depicting the levels of MET in Ecad high and Ecad low subsets of MBM and EM (left panel, p = 1.4e-04/p = 0.41) or in all subtypes of MBM and EM (p = 2.7e-05) in high and low proliferating tumor cell subsets (right panel, p = 9.1e-03).d Comparative principal component (PCA) representation of primary tumors (PT), extracranial metastases (EM) and MBM (MBM_TCGA) of the TCGA-SKCM cohort as well as MBM of our study (MBM_CHA, EGAS00001005976) depicting gradual levels of MET and Ecad (CDH1) expression.The panels below show a comparison of levels of MET, Ecad, MITF and NGFR of selected tumors showing distinct and overlapping cell states.Expression levels (log2, FPKM) are color coded.e IHC of selected MBM for MET and MITF validated the two subsets.f, g Expression and activation status of MET in BRAF wildtype (wt; Pts 14, 36) and BRAF V600E/R mutated MBM (Pts 28, 29, 31).Phosphorylation of MET at residues Y1234/1235 is critical for kinase activation.h IHC of indicated tumors for co-localization of pMET Y1234/1235 (brown) and Iba1 (red) demonstrating potential activation of MET receptor signaling tumor cells by stromal cell-secreted HGF.i Heat map representing expression levels of regulators and targets of interferon signaling and immune-related genes showing clustering according to the level of ITGB7 expression.Box and whisker plots show median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (c, d).Significance was determined by unpaired, two-sided t-test (c, d) V600K , ~50 % p.BRAF N581Y , 48 % / p.NRAS G12C , 69 , IC 50 =2.3µM BMC1-M4, IC 50 =1.7 µM BMC2, IC 50 =2.2µM A375, IC 50 =1.4µM MeWo, IC 50 =2.5 µM Dose [nM]

Fig. 5
Fig.5 Inhibitors of MET receptor decrease the growth of brain metastatic and conventional melanoma cell lines.a Comparative IHC of selected MBM for levels of phosphorylated and activated MET receptor (pMETY1234/1235 ) and ribosomal protein S6 (pS6 235/236 ) of consecutive sections suggesting MET-associated activation of mTOR signaling.b Immunofluorescence microscopy of lymph node-metastatic (T2002) and brain metastatic (BMC53) patient-derived melanoma cell lines for the co-occurrence of MET (red) and pS6 235/236 (green).DAPI served as a nuclear counterstain.c qPCR analysis of BMCs for expression of MET receptor, bars indicate median levels ± SD of three biological replicates.d Gross initial ARQ197 sensitivity test of BMC53 and BMC1-M1 cells showing high and low levels of MET expression.Cell density was determined by crystal violet staining.e Broad range determination of sensitivity of BMCs, T2002 and conventional melanoma cell lines (A375, A2058, MeWo) to METi PHA-665752 and ARQ197.Cell density and BRAF mutation status are indicated.Dotted line depicts the estimated range of IC 50 .f PHA-665752 dose-response fit curve-based calculation of IC 50 values of A375 cells with overexpression of NGFR or RFP control cells and MeWo cells.g Dabrafenib dose-response fit curve-based calculation of IC 50 values of BMCs exhibiting different BRAF mutations (BMC2 p.N581Y , BMC4 p.V600K ) and A375 p.V600E , A2058 p.V600E cells.h, i Live cell imaging-based tracking of confluence (%) of BMC2 and BMC4 cells in dependence of increasing doses of ARQ197.Shown are median values ± SD of eight technical replicates.A representative out of two experiments is shown.j ARQ197 dose-response fit curves of BMCs, T2002 and conventional cell lines.Calculated IC 50 values are indicated, suggesting the response of dabrafenib-resistant cell lines to METi.k Bar diagram summarizing IC 50 values (nM) indicating the response of indicated cell lines to ARQ197.The BRAF status is color coded.l Working model suggesting the activation of MET receptor signaling in adjacent tumor cells and in (reactive) microglia (RM) by microglia-released HGF.IRF-mediated HGF expression in turn is triggered by immune cell (monocytes/macrophages, M) released interferon-gamma.HGF binding to tumor cell (TC) expressed MET receptor directs downstream activation of the RAS/ RAF/MEK/ERK and the PI3K/AKT/mTOR/p70S6K branch.The latter is leading to phosphorylation and activation of the ribosomal protein S6.Box and whisker plots show the median (center line), the upper and lower quartiles (the box), and the range of the data (the whiskers), including outliers (c)

Fig. 6
Fig. 6 HGF/MET receptor signaling might be activated in tumor cells at immune cell/TAM dense areas.Schematic representation of our working model suggesting the interaction of tumor cells with stromal cells, particularly microglia and immune cells, consequentially lead-